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Abstract Linear Elastic Fracture Mechanics (LEFM) 
provides a consistent framework to evaluate quantita- 
tively the energy flux released to the tip of a growing 
crack. Still, the way in which the crack selects its veloc- 
ity in response to this energy flux remains far from com- 
pletely understood. To uncover the underlying mecha- 
nisms, we experimentally studied damage and dissipa- 
tion processes that develop during the dynamic failure 
of polymethylmethacrylate (PMMA), classically con- 
sidered as the archetype of brittle amorphous materi- 
als. We evidenced a well-defined critical velocity along 
which failure switches from nominally-brittle to quasi- 
brittle, where crack propagation goes hand in hand with 
the nucleation and growth of microcracks. Via post- 
mortem analysis of the fracture surfaces, we were able 
to reconstruct the complete spatiotemporal microcrack- 
ing dynamics with micrometer/nanosecond resolution. 
We demonstrated that the true local propagation speed 
of individual crack fronts is limited to a fairly low value, 
which can be much smaller than the apparent speed 
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measured at the continuum-level scale. By coalescing 
with the main front, microcracks boost the macroscale 
velocity through an acceleration factor of geometrical 
origin. We discuss the key role of damage-related in- 
ternal variables in the selection of macroscale fracture 
dynamics. 

Keywords Dynamic fracture • amorphous polymers • 
fracture energy • fractography • high-deformation rate • 
microcracks ■ PMMA 



1 Introduction 

Dynamic crack propagation drives catastrophic mate- 
rial failure and is usually described using the Linear 



Elastic Fracture Mechanics (LEFM) framework (Fre- 



und 1990 Ravi-Chandar 2004). This theory considers 



the straight propagation of a single smooth crack and 
assumes all energy dissipating processes to be entirely 
localized in a small zone at the crack tip, so called frac- 
ture process zone (FPZ). Crack growth velocity v is 
then selected by the balance between the mechanical 
energy that flows within the FPZ per time unit and the 
dissipated energy within the FPZ over the same time 



unit. This yields (Freund 19901: 



r~(l-v/C R )K 2 (c)/E, 



(1) 



where Cr and E are the Rayleigh wave speed and the 
Young modulus of the material, respectively, r is the 
fracture energy, and K(c) is the Stress Intensity Factor 
(SIF) for a quasi-static crack of length c. K only de- 
pends on the applied loading and the cracked specimen 
geometry, and entirely characterizes the stress field in 
the vicinity of the crack front. 
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LEFM predictions agree well with observations as 



long as the crack growth is sufficiently slow ( Bergkvist 



1974 1. However, large discrepancies are reported at high 



speed (Ravi-Chandar 2004 Bouchbinder et al 20101. 



In particular, the maximal crack speeds attained exper- 
imentally in amorphous materials are far smaller (typ- 
ically by a factor of two) than the limiting speed Cr 
predicted by Eq. [T] 




Imposed 
Displacement rate 



Support Rollers 



The existence of a micro-branching instability (Finet erg 



et al 1992 I at a critical velocity Vb (typically of the or- 
der of 0.35 — 0.4(7r) permits to explain part of this 
discrepancy: Beyond u&, the crack front splits into a 
multiple crack state that cannot be described by Eq. IT 
anymore. However, it is commonly stated (Sharon and 



Fineberg 1999) that LEFM works well below «&, or 



more generally when microbranching instabilities are 
absent. As a matter of fact, most recent works focused 
on failure regimes beyond Vb (see e.g. ( jGumbsch et al[ 
19971 |Adda-Bedia~et~aIl [19991 jHenry and Levinc | 2004 



Bouchbinder et al 2005 Spatschek et al 96 Henry 



2008 ) ) . Still, several observations reported at lower speeds 
remain puzzling. In particular, even for velocities much 
lower than Vb, the measured dynamic fracture energy is 
generally found to be much higher than that at crack 



initiation (Sharon and Fineberg 1999 Kalthoff et al 



1976 Rosakis et al 1984 Fond and Schirrer 2001). 



The experiments reported here were designed to bet- 
ter understand the mechanisms that select crack veloc- 
ity in dynamic fracture. In this context, we developed 
an experimental setup (presented in section 2) that per- 
mits to characterize over a wide range of crack speeds, 
in particular speeds below Vb, the dissipative and dam- 
age processes that develop in polymethylmethacrylate 
(PMMA), classically considered as the archetype of brit- 
tle amorphous materials. This setup allowed us to evi- 
dence a novel critical velocity, smaller than Vb, beyond 
which crack propagation goes hand in hand with the nu- 
cleation and growth of microcracks ahead of the main 
crack front (section 3) . Via accurate post-mortem anal- 
ysis of the patterns let on fracture surfaces, we were 
able to reconstruct the full spatio-temporal dynamics 
of these microfailure events (section 4). Their statis- 
tics and their dependency with crack tip loading have 
been characterized (section 5). In section 6, we will show 
how microcracking acts and selects the apparent veloc- 
ity measured at the continuum-level scale. 

2 Experimental setup 

Fracture tests were performed in the wedge splitting 



geometry (Bruhwiler and Wittmann 1990 Karihaloo 



and Xiao 2001 1 depicted in Fig. 1. Specimens were pre- 




pared from PMMA parallelepipeds of size 140 x 125 x 



Fig. 1 Sketch (top) and photo (bottom) of the experimental 
setup. 



15 mm 3 in the propagation (ar-axis), loading (y-axis), 
and thickness (z-axis) directions, respectively. Table 1 
gives the main characteristics of the PMMA used in 
this study. Subsequently, a notch was shaped (i) by cut- 
ting a 25 x 25 mm 2 rectangle from the middle of one 
of the 125 x 15 mm 2 edges; and (ii) by subsequently 
adding a 10-mm groove deeper into the specimen. Two 
steel blocks equipped with rollers were then placed on 
both sides of this notch and the specimen was loaded by 
pushing a wedge (semi-angle of 15°) at a small constant 
velocity (40 /xm/s) between these two blocks. This per- 
mits (i) to spread the loading force over a large contact 
area and prevent any plastic deformation of PMMA at 
the loading contacts; and (ii) to suppress friction in the 
system. As a result, the vicinity of the crack tip can be 
assumed to be the sole dissipation source for mechanical 
energy in the system. 

In such a wedge splitting geometry, the SIF de- 
creases with crack length. To obtain dynamic failure, 
we then introduced a circular hole of tunable diameter 
from 2 mm to 8 mm at the tip of the seed crack. This 
hole blunts the seed crack and delays the propagation 
onset. This increases the mechanical energy stored in 
the sample at the onset of crack propagation. It allowed 
us to widen the range of SIF (from 0.9 MPa.-yrn to 
4.5MPa.y / m) and that of crack velocity (from 75 m/s s=s 
0.08Cr to 500 m/s sw 0.57Cr) accessible in the experi- 
ments. 

The apparatus used to carry out the experiments 
is homemade (Fig. [TJ) . One side of the machine is fixed 
while the other moves via a stepper motor (Oriental mo- 
tor EMP400 Series) allowing incremental displacements 
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Fig. 2 Measurement of instantaneous crack velocity, a: Voltage as a function of time as recorded via the oscilloscope, b: Zoom 
that permits to distinguish the individual jumps yielded by the successive cuts of the conductive lines as crack advances, c: 
Resulting position of crack tip as a function of time, d: Final curve showing the variations of crack velocity with crack length. 
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Fig. 3 Computation of SIF. Top: Typical mesh used for finite elements calculations, in order to access the stress/strain 
fields in the experiments. Red: polymeric layer. Black: sample. Blue: L-shaped block. Green line: cracked line. The slide link 
connected to the L-shaped block is used to model the motion of the contact point between the pushed wedge and the roller. 
Specimen loading is achieved by translating horizontally the slide link, from (0) to (1), over a distance u we d ge selected so that 
the horizontal force applying on the slide link is half that measured experimentally. Bottom-Left: Zoom on the meshing in the 
transition region (red) between coarse meshing (1 mm mesh size) in the bulk and fine meshing close to the crack tip. On the 
right is shown part of the circular hole at the seed crack tip. Bottom-Right: Zoom on the crack-tip region (green), meshed 
with a size of 1 nm. 
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£(GPa) 



Density (kg/dm 3 ) Co, (m/s) C s (m/s) Cr (m/s) 



PMMA 2.8 ±0.2 0.36 ±0.08 



1.18 ±0.02 



2010 ±60 950 ±30 



)±20 



Table 1 Mechanical properties of the PMMA used in our fracture experiments. E andf denote the Young modulus and the 
Poisson ratio, respectively. Cd , Cs , and Cr denote the speeds of dilational, shear and Rayleigh waves respectively. 



as small as 40 nm. The compressive force is measured, 
up to 20 kN via a S-type Vishay load cell 363 Series. 
In all performed tests, both the displacement and the 
force are recorded via a computer with a time resolu- 
tion of 1 s. In particular, the value of the applied force 
at propagation onset is recorded and later used in the 
finite element analyses described thereafter. 

To measure the instantaneous crack velocity, we used 
a modified version of the electrical resistance grid tech- 



nique (Cotterell 



1968 



Anthony et al 



et al 1996 Fineberg et al 1992 ) . A series of 90 paral 



1970 Boudet 



lei conductive lines (2.4-nm-thick Cr layer covered with 
23-nm-thick Au layer) were deposited on one of the two 
140 x 125 mm 2 sides of the specimen using magnetron 
sputtering deposition through a polymer mask glued 
on the sample surface. Stripes width and length are 
500/xm and 100 mm, respectively, and they are sepa- 
rated by gaps of 500/xm. Once the sample is installed 
in the compression machine, each line is connected in 
series with a resistor of 30 ki? - this value was chosen 
to be extremely high with respect to that of the strip 
(about few ohms) so that the latter is negligible. As 
the crack propagates throughout the specimen, it suc- 
cessively cuts the conductive lines, which, each time, 
makes the global system resistance increase by a con- 
stant increment. The time locations of these jumps are 
detected via an oscilloscope (acquisition rate: 10 MHz) 
through a voltage divider circuitry. In order to increase 
the accuracy of jump detection, we used the four chan- 
nels of the oscilloscope in cascade with different offsets, 
so that each channel is triggered as the previous one 
gets out of its range (Fig. [2] A and B). As a result, we 
obtained the crack length c(t) as a function of time with 
space and time accuracies of 40 /im and 0.1/Lts, respec- 
tively (Fig. J2JC). The profile of instantaneous velocity 
v(c) as a function of crack length can then be easily 
deduced (Fig.|2JD). 

Finite element analysis was used to estimate the 
SIF evolution during failure. Figure |3]Top represents 
the meshing of the complete system. The average mesh 
size is 1 mm and reduces logarithmically down to 1 nm 
at the crack tip (Fig. |3JBottom), in order to resolve in 
space the crack tip opening. Only half of the system 
is required for symmetry reasons. Young modulus and 
Poisson ratio are set to 2.8 GPa and 0.36, respectively 
(see Tab. 1). The boundary conditions are the follow- 
ing: (i) The left edge is perfectly adhesive to a polymeric 



layer of thickness 5 mm, Young's modulus 3 GPa and 
Poisson ratio 0.41; (ii) The left edge of the layer has a 
no-displacement condition; (iii) On the right, the notch 
is loaded through a L-shaped block, of thickness 5 mm, 
Young modulus 300 GPa and Poisson ratio 0.4, per- 
fectly adhering to the sample; (iv) The L-shaped block 
is bounded to a virtual roller, the center of which can 
only move along a line parallel to the wedge side, at a 
distance equal to the roller radius; (v) The crack edge 
is stress free; and (vi) Vertical displacement is forbid- 
den along the uncracked part of the symmetry plane 
(top line of Fig. |3JTop). For each sample, the equilib- 
rium position of the wedge yielding the measured ap- 
plied load at the onset of crack propagation was de- 
termined using a plane stress static finite element code 
(Cast3M). Quasi-static crack propagation is then sim- 
ulated by increasing the crack length while imposing 
the wedge position to remain constant. This latter as- 
sumption is experimentally justified. Indeed, it has been 
observed that the load cell signal started to be modi- 
fied only after a few hundreds of microseconds (typical 
time of an experiment) after crack initiatiorF] This is 
attributed to the time required for the sound waves to 
travel from the crack tip to the load cell. As a result, 
the wedge position can be assumed to be practically 
constant during crack propagation. Static SIF is then 
determined, for each value c of crack length, using the 



J-integral method (Rice 1968). The validity of our SIF 



calculation was benchmarked against literature results 
for wedge-splitting configurations: Our code agreed to 
better than 2% with the results of IKarihaloo and Xiaol 



(20011 when using their system parameters. 



To image the post-mortem fracture surfaces, we used 
a Leica DM2500 microscope. Most of the images were 
taken with a x5 objective under polarized light. The 
resulting images are 1.4 x 1 mm 2 in area, with a pixel 
size of 677 nm. We also imaged the fracture surfaces 
with an interferometry profilometer (FOGALE Nan- 
otech) which allows one to gather topographical infor- 
mation. These 3D images were taken with a x5 ob- 
jective. Their area is 1.4 x 1 mm 2 with a pixel size of 
1.86 /im. 



1 To measure this time shift, we directly connected the 
Wheatstone bridge of the load cell to an oscilloscope, without 
passing through the signal conditioner. This latter, indeed, 
imposes a time resolution of 1 s. 
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Fig. 4 Variation of SIF K (left), instantaneous crack velocity v (center) and fracture energy r (right) as a function of crack 
length c for five experiments with different stored mechanical energy Uq at propagation onset: 2.0 J (A, A' and A"), 2.6 J (B, 
B' and B"), 2.9 J (C, C and C"), 3.8 J (D, D' and D"), and 4.2 J (E, E' and E"). 
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Fig. 5 Fracture energy f as a function of crack velocity v 
for five different experiments with different stored mechanical 
energies Uq at crack initiation. The two vertical dashed lines 
correspond to v a and vi,. The two horizontal dashed lines in- 
dicate the confidence interval (95%) for the me asured fracture 
energy K^/E at crack initiation (Taken from |Scheibert et al| 



(2010)) 



3 Selection of fracture energy with velocity and 
evidence of a brittle/quasibrittle transition 

In each of the experiments performed, both the SIF 
profile K(c) and the instantaneous velocity profile v(c) 
have been determined (Figs. EJLeft and center). From 
these curves, one can derive the profiles of fracture en- 
ergy, r(c), using Eq. [I] (Fig. EJRight). We note here 
that Eq. [T]is strictly valid in the case of crack propaga- 
tion in an elastic half-space, i.e. in the absence of waves 
emitted by the crack and coming back to it after reflec- 



tion on sample boundaries (Freund (19901; Goldman 



et al (2010 1). We argue that our experiments are es- 
sentially unaffected by such reflected wave for at least 



two reasons: i) as measured by Boudet et al (19951, 
the sound energy emitted by a crack in PMMA lies 
within the range 1 — 4 J/m , a very small value com- 
pared to the energy required to fracture the material 
(a few kJ/m , see below); and ii) this radiated energy 
is quickly dissipated within the material's bulk - The 
acoustic attenuation coefficient has been measured to 
be 0.67 dB/(cm.MHz). 

Figure [5] superimposes the various curves r vs. v 
measured in the various experiments performed at var- 
ious levels Uq of the initially stored mechanical energy. 
Several regimes can be evidenced: 

— For small v, r roughly remains constant, close to 
K^/E, as expected within standard LEFM; 

— As v increases and reaches a first critical velocity 
v a « 165 m/s ~ 0.19Cr, r suddenly increases to 
a value about four times larger than r(v = 0) = 
Kl/E- 

— Beyond v ai r slowly increases with v; 



— As v reaches a second critical velocity Vb ~ 350 m/s ~ 
0.36(7r, r starts to increase rapidly again. It seems 
to diverge at v « 450 m/s ps 0.5Cr. 

Note also the fairly good collapse of the different 
curves r(v) below vj, and the large dispersion above. 
This suggests that Eq. [T] is relevant for v < Vj,, - pro- 
vided a suitable velocity dependence r(v) is ascribed-, 
but not beyond. The second critical value Vb ~ 350 m/s R 
0.36(7r is found to correspond to the onset of micro- 
branching instability widely discussed in the literature 



(see e.g. (Fineberg et al 1992 Sharon and Fineberg 



1999)). The first critical value v a « 165 m/s ps Q.19C R 
was observed in our series of experiments and reported 



in (Scheibert et al 2010) for the first time. This observa- 



tion was made possible by the use of the wedge-splitting 
geometry, which is a decelerating crack configuration. 
It offers many data points at relatively low velocities, 
contrary to most other devices used before. The rapid 
increase of r with v around v a provides a direct inter- 
pretation for the repeated observations of cracks that 
span a large range of r, while keeping a nearly con- 



stant velocity about 0.2C^ (see e.g. ( Ravi-Chandar and 



Knauss 1984 Ravi-Chandar and Yang 1997 1). 

To shed light on the nature of the transition at 
v = v a evidenced on the curve relating f to f, we exam- 
ined the post-mortem fracture surfaces and their evolu- 
tion as crack speed increases (Fig.[6l. For v smaller than 
v a , fracture surfaces remain smooth at optical scales. 
But above this threshold, conic marks start to be ob- 
served. These marks remain confined within a rough- 
ness scale of few micrometers at the fracture surfaces, 
contrary to the microbranches that deeply develop into 
the specimen bulk above v\> (see side views and topo- 
graphical images in Fig. [71. 

Similar conic marks were reported in the fracture of 



other brittle materials, among which polystyrene (Regel 



19511, oxide glasses (Smekal 1953 Holloway 1968), cel- 



lulose acetate (Kies and abd G. R. Irwin[ 19501, poly- 



crystalline materials (Irwin and Kies 1952), or homalitc 



associated ( Smekal 



Rabinovitch et al 



(Ravi-Chandar and Yang 19971. They are classically 



1953; Ravi-Chandar and Yang 1997 



2000| with the presence of microc- 



racks that nucleate and grow ahead of the main crack 
front, and subsequently coalesce with it. 

Figure [8] shows the surface density p of conic marks 
as a function of crack speed v. Below v a , no mark is ob- 
served, irrespectively of the chosen magnification (from 
x5 to x50). Above this value, p increases almost lin- 
early with v — v a . The precise correspondence between 
the critical velocity v a at which the curve r(v) ex- 
hibits a kink and that at which conic marks start to 
be observed suggests that both phenomena are the sig- 
nature of the same transition. Above v a , PMMA failure 
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Fig. 6 1 x 1.4 microscope images of fracture surfaces created at various speed v. 
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Fig. 7 Profilometric 1.4 x 1 mm 2 surface (Top) and optical 0.7 x 0.25 mm 2 side-view (Bottom) images taken at velocity 
smaller than v a (left), between v a and uj, (center), and above vt (right). Crack has propagated from left to right. 



switches from nominally brittle to quasi-brittle and goes 
hand in hand with microcracking that develop ahead of 
the main front. In the following, we will use the term 
"Damage Zone" (DZ) to refer to the zone where micro- 
cracks develop in the vicinity of the main crack tip. We 
will distinguish this zone from the FPZ (smaller than 
the DZ) that embeds dissipative mechanisms (crazing 
for instance) at the tip of each (micro)crack front. 



Note that the region of the fracture surface that 
bears conic markings also has a characteristic aspect to 
the naked eye, as it scatters ligth very efficiently. This 
is markedly different from the high reflectivity of the 
region that bears no microcrack. In fractography these 
regions are classically referred to as the mist and mirror 
zones, respectively (see e.g 



(Hull 



1999 Rabinovitch 



and Bahat 2008| )). Our observation therefore strongly 

suggests that the mirror/mist transition is simply the 
morphological counterpart of the brittle/quasi-brittle 
transition that occurs at v a 
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Fig. 8 Surface density p of conic marks as a function of 
crack velocity for five different experiments with different 
stored mechanical energies Uo at crack initiation (Taken from 
|Scheibert et al| ( |2010] >) . 
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Fig. 9 A: Typical microscope image of the post-mortem fracture surface. Bright regions indicate the nucleation centers. B: 
Interpretation of the conic marks: Red dots sketch the successive positions of the fronts that, by interacting, will be giving 
rise to the green branch of the conic mark. These two microcracks form at t = and t = T, then radially grow at speeds C\ 
and C2- Branch fitting via Eq.plpermits to measure C2/C1. Note that a conic mark can be made of several branches. C: When 
c-2 = C\ , branches are (mathematical) true conies and the nucleation distance d n between the front of the incident microcrack 
(1) and the nucleation point of the forming one (2) is twice the apex-focus distance O2P12 (Taken from|Guerra et al (2012)). 



4 Practographic reconstruction of individual 
microfailure events 



In order to further characterize and understand the dy- 
namics of microcracking events which develop during 
fast fracture in PMMA (v > v a ), we proposed a numer- 
ical reconstruction based on the post-mortem analysis 



of the traces left by these events (Guerra et al 2012 1. 



The typical time between two successive microcracking 
events is about 10 ns, which makes them inaccessible us- 
ing standard techniques as e.g. fast imaging or acoustic 
emission analysis. 

The first step is to identify where microcracks have 
initiated. On the optical images of post-mortem frac- 
ture surfaces (Figs.[9]A and B) , many well defined points 
of high optical reflectivity can be seen. These bright/white 
areas are believed to come from large plastic defor- 
mations accompanying the nucleation of microcracks. 
Thus, they allow us to determine precisely the position 
of the different nucleation sites. For some microcracks, 
the precise location of their nucleation site is further 
constrained by the convergence of fragmentation lines 
onto it. 



The second step is to infer, from the conic pat- 
terns, the velocities at which the various microcracks 
have grown. It is commonly admitted in the litera- 



ture) Smekal 1953 Ravi-Chandar and Yang 1997 Ra 



binovitch et al 20001 that the conical marks indicate 



the intersection points between two interacting microc- 
racks. Let us then consider two microcracks that prop- 
agate radially along slightly different planes with veloc- 
ities c\ and C2 (Fig. [9JB). The intersection between the 
initial (micro)crack front and the secondary microcrack 
leaves a visible trace (i.e. a tiny height difference) on 
the fracture surfaces. In the coordinate system (e x ,e a ) 
centered on the site of nucleation of the first microcrack 
and chosen so that the x-axis passes through the two 
centers of nucleation, the equation describing the mark 



left postmortem can be written as (Guerra et al 2012 1: 



= ± 



2c 2 CiT 


(c 2 -l) 2 A 


-d) 2 + ( 


c* + l 



'(,_ D ( *_ 1)+ „(-)■■ 



(c 2 



(c 2 - l) 2 



m 



2x 

~A 



-1I/2 



(2) 
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Fig. 11 Reconstructed sequence of crack propagation and 
microfailure events at the microscale for two different val- 
ues of K, namely K = 2.53 MPa.ymi (left) and K = 
4.18 MPa.ymi (right). The blue part corresponds to un- 
cracked material and we ascribed an arbitrary gray level to 
each of the growing microcrack. From this, one can deduce the 
intersection points between microcracks that coincide with 
the observed conies marks on the fracture surfaces. The size 
of the two zones of analysis is 2.5 x 2.5 mm 2 . The macroscopic 
crack propagates from left to right. 



where c = C2/C1, A is the distance between the two nu- 
cleation sites, and r is the time interval between two mi- 
crocracking events. The analysis of this equation shows 
that the number c\T / A sets the aspect ratio of the frac- 
tographic trace, while the ratio c = c^/ci fixes its form 



?*»^ ''Vi»1 ^-r-' 


mm ji . .^^& 


r-^A 




s 300 ym i 








Fig. 12 Qualification of the procedure: Superimposition of 
the conic pattern as obtained from reconstruction onto that 
observed by fractography. A: K = 2.77 MPa.y'E. B: K = 
4.18 MPa. v / m. Red dots indicate nucleation cen ters (T aken 
from the supporting information of|Guerra et al| (|2012[|). 



We used Eq.[2]to directly extract the ratio c = c 2 /ci 
for each couple of interacting microcracks (~ 400 cou- 
ples per image) on different optical pictures obtained at 



different speeds (Fig. 10 1 . Regardless of the macroscopic 
crack velocity, we obtain for c a Gaussian distribution 
centered on 1 with a standard deviation of about 0.03 



(see Guerra et al (20121 for details) 



(Fig. 10 1. This analysis permits to demonstrate that all 
microcracks grow at the same velocity c m inside the 
DZ. Note that, up to now, nothing prevents c m to vary 
with v or K . The mechanism that fixes the value for c m 
will be discussed later in this paper (see section pi . By 
using C2 = c\ — c m in Eq. [2] the following relation can 
be obtained: 
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Fig. 10 Direct extraction of the relative speed between two interacting microcracks at three different microcrack densities. 
(Top) Zones of investigation. Each conic branch has been attributed a given color and the nuclei of the two corresponding 
interacting microcracks has been joined by a dotted segment of the same color. Note that a conic-like mark is often made of 
several of these conic branches. The ratio C2/C1 is the only adjustable parameter in equation Eq. [2] to determine the branch 
geometry once the nuclei position and the branch apex are set. (Bottom) Corresponding distributions for C2/C1.: In the three 
cases, the distributions are found to fit normal distributions of mean value 0.98 — 0.99 and standard deviation 0.03 — 0.04, 
irrespective of p (Taken from the supporting information of|Guerra et al|||2012[)). 



A 



= ± 



d n (2A - d n ) 

{A-2d n y 

d n (2A - d n ) 

(a- 2d n y 



x\ 2 d n (2A-d n ) / x^ 
AJ (A-2d n ) 2 \A 



-4 



-1 1/2 



(3) 



where d n = (A — c m T) represents the critical distance 
between the front of the incident microcrack and the 
center of the nucleated microcrack at the instant of its 
nucleation. This distance is twice the distance between 
the apex and the focus of the considered microcrack 
and, as such, d n can be directly extracted from fracto- 
graphic images (Fig. [9JC). 

We are now in a position to reconstruct the full dy- 
namics of microcracking events from the knowledge of 
two series of parameters that can be directly extracted 
from fractography: The position of nucleation sites and 
the critical nucleation distance d n associated to each 
conic mark. We analyzed several areas of the fracture 
surfaces, corresponding to various macroscopic veloc- 
ities above the microcracking onset (v > v a ). In all 
cases, nine optical images with partial overlap were 
recorded and gathered into a single large image - This 
ensures adequate statistics in the following analyses. 
The methodology to extract the data is the following: 

— We record the coordinates x and y of the nucleation 
sites; 



— We estimate the nucleation distance, d n , by infer- 
ring, for each nucleation site, who was the parent 
microcrack from the relative position of the conies 
apex with respect to the nucleation site (family cri- 
terion). Then, d n is twice the apex-nucleation center 
distance. 



Following Ravi-Chandar and Yang (1997), the re- 



construction is initiated with an originally straight crack 
front positioned at the left of the zone of analysis. This 
front propagates from left to right, at a constant ve- 
locity of 1 pixel/timestep. Once the distance between 
one of the microcrack nucleation sites and the primary 
crack front reaches d n , a secondary microcrack starts to 
grow radially at the same velocity. The front is now the 
combination of both the translating straight front and 
the radially growing circular front. Beyond a given time, 
these two fronts intersect and define the position of the 
conic mark. When the growing front reaches again a 
distance d n associated with another site, a new micro- 
crack is nucleated. The procedure is repeated until all 
the nucleation sites have been triggered. 

Figure [TT] shows two examples of reconstruction for 
two different macroscopic velocities and Fig. [12] super- 
imposes the conic marks obtained by the simulation to 
those really observed in fractography. We note the fairly 
good agreement between experiment and reconstruc- 
tion away from the edges (in the center of the images) . 
This permits to validate the method. Thanks to it, we 
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are able to reconstruct deterministically the dynamics 
of micro-damage during fast fracture in PMMA from 
the observation and analysis of the patterns observed 
on the post-mortem fracture surfaces. We emphasize 
the resulting resolution of ~ 1 /im (1 pixel) in space, 
and of 10 ns in time (pixel size/c m with c m ~ 200 m/s 
as it will be demonstrated further in the paper). This is 
far better than what can be obtained via conventional 
experimental mechanics methods as e.g. acoustic emis- 
sion or fast imaging methods. 

Note however some discrepancies between experi- 
ments and reconstructions on the left, top and bot- 
tom edges of the analyzed zones. Such finite size prob- 
lems are unavoidable. Along the top and bottom sides 
of the image, the influence of microcracks outside the 
field of view are naturally ignored. Also, a non-realistic 
straight vertical front has been used, as we cannot pre- 
dict the precise instants at which the leftmost centers 
have turned to microcracks during the reconstruction. 
In order to test the sensitivity of the reconstruction to 
the initial front shape, we have performed different sim- 
ulations with the same inputs except the initial front 
shape. In Fig. 13 three different cases are tested : (i) a 



straight vertical front, (ii) a sinusoidal front with an am- 
plitude and a period equivalent to the mean value of a 
wavy front far from the first step in previous reconstruc- 
tion and (iii) the same sinusoidal front vertically shifted 
by half a period. We can observe that, except of course 
at the very beginning of the simulation (i.e. on left side 
of the images) , the difference between these three cases 
are very small. Moreover, these small differences tend 
to disappear as the front to propagate, i.e. as more and 
more microcracks are involved. In some cases, we have 
observed differences propagating over the whole image. 
Perfect reconstructions would be reached only by ana- 
lyzing the whole fracture surfaces, and not some partial 
areas. Knowing all the history of crack propagation is 
the only solution for the determination of the real initial 
front shape which could allow a perfect reconstruction. 



Fig. 13 Influence of initial crack front morphology on the re- 
construction for p = 45.0 mm -2 (K = 3.65 MPa.-y/m). Three 
different initial conditions were used. (A) Straight vertical 
line. (B) Vertical sinusoidal shape with a period of 186 /im 
and a peak-to-peak amplitude of 242 //m. (C) Same sinusoidal 
shape, but translated vertically over half a period (Taken from 
the supplementary information of|Guerra et al||2012)). 



5 Quantitative selection of the microdamaging 
state by crack loading 

The reconstruction method, described in the previous 
section, shows that the dynamics of micro-damage is 
entirely determined by the position of nucleation sites, 
the nucleation distances d n , and the microscopic veloc- 
ity of the fronts c m . To understand what selects these 
three parameters, they will be explored in more detail. 
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Fig. 14 Cumulative distribution of the number of nucleation 
centers contained in square regions of lateral size a, for two 
different fractographic images. Solid lines represent Poisson 
fits of parameter pa 2 , where the fitted value p can be identified 
with the center density. 




Stress intensity factor (MPa.m ) 



Fig. 15 Variation of the density p of nucleation centers with 
applied SIF. Black disks correspond to the eight images on 
which post-mortem reconstruction was performed. Red line is 
a fit via Eq. R] K a and Kt are associated with microcracking 
and microbranching onset, respectively (Taken from |Guerra| 
|et~al|p0l2) ). 



5.1 Spatial distribution of nucleation centers 

To characterize the spatial distribution of nucleation 
sites, we computed the distribution of the number of 
sites that can be observed on a square of a given size 
a. Figure [14] shows these distributions for two different 
fractographic pictures, i.e. for two values of K. These 
distributions follow a Poisson distribution with a single 
parameter pa 2 where p is constant over a given fracto- 
graphic picture but varies with K . This reflects a ho- 
mogeneous and uncorrelated distribution of the sites all 
along the post-mortem fracture surface. Then, p can be 
identified with the surface density of nucleation sites 
and fully characterizes their spatial distribution. 

Figure [15] presents the evolution of p with loading 
K . The density is naturally zero for K smaller than K a 
(SIF value at microcracking onset, i.e. when v = v a ) 
as there is no microcracking. Beyond this value, it in- 
creases with K. Also, this curve exhibits large fluctua- 
tions and seems to saturate when the density reaches a 
value p sat ~ 45.5 mm 2 , beyond K^. Note that this latter 
value coincide with microbranching onset and, therefore 
the p(K) curve does not really make sense above K\, 
since the front has splitted into various microbranches, 
which prevents us from correlating the macroscopic K 
value to the number of conic marks on the main crack 
only. 



To account for the curve p(K), let us consider ( Scheib- 



ert et al 2010 1 that a discrete population of weak lo- 



calized zones is present within the material with a bulk 
density p v . Let us then assume that a microcrack can 
nucleate in these areas provided the two following con- 
ditions are fulfilled: 

— The local stress at the considered zone is large enough, 
i.e. larger than a given threshold value a* (smaller 
than the yield limit ay of the material) ; 

— The considered zone is far enough from the main 
crack front to allow the microcrack to grow, i.e. at 
a distance greater than d a . 

The density p of nucleated microcracks per unit area 
of fracture is then given by p v {h± — d a } where h± is 
the thickness of the layer where the stress has exceeded 
a* (h± is measured perpendicular to the mean plane of 
fracture) . The singular form of the stress fiel d at crack 



tip leads to write h± ex K /a* . which yields (Scheibert 



et al 20101: 


for 
for 




p = 
p ex K 2 - Kl 


K<K a 
K>K a 



(4) 

with K a tx a* \fda. This equation, plotted in red on Fig. 
|15| reproduced quite well the experimental data below 
K b . 
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Fig. 16 Cumulative distribution of the nucleation distances 
d n , as measured for two different fractographic images. Solid 



lines are fits with P(d n ) = ((d„ 



d n )/(dr, 



d r , 



O)* 



As the determination of p from fracture surface is 
much more accurate and consensual than the determi- 
nation of K (based on FEM calculation with inherent 
approximations), in the following, we will use p rather 
than K as the control parameter as a function of which 
the various quantities further defined in this paper will 
be plotted. Indeed, p fully characterizes the spatial dis- 
tribution of nucleation sites and monotonically evolves 
with K within the relevant range K a < K < Kb- 



5.2 Chronology of microfailure events 



Figure 16 shows the cumulative distribution of nucle- 
ation distances d n measured on two typical fracture sur- 
faces. The probability density function (derived from 
the cumulative distribution) is zero for d n < d m i n , max- 
imum at d min , subsequently linearly decreases, and is 
finally zero as d n exceeds d max . This distribution form 
was observed for all the fracture surfaces, irrespectively 
of p. d max is found to increase with p. while d m ; m de- 
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Fig. 17 Top: Mean distance at nucleation d n as a function of 
p. the red curve is a linear fit over p < p sa t = 45.5 mm and 
a plateau at d n sa = 50 pm over p 3a t- Bottom: Saturation 
of d n and avalanches. The red line is the same as in the top 
figure, ft is compared to the mean nearest-neighbor distance 
in a Poissonian distribution (black solid line). Solid squares 
indicate the proportion of microcracks involved in avalanches, 
as computed from the reconstruction. 



creases with p and remains very close to zero (Guerra 



et al 2012). 



The variation of the mean value d n of d n with p 
exhibits two regimes (Fig. [l7JTop): (i) A linear increase 
with p followed by (ii) a saturation (d n ~ 50 pm) 
when p is greater than the value p sat associated in Fig. 
15 with the microbranching onset (at K = Kb). 

To understand the saturation origin, we plotted, 
on Fig. [l7]Bottom, the evolution of the mean nearest- 
neighbor distance (Ar) with p as it is predicted for a 2D 
Poissonian spatial distribution ((Ar) = 1/2^/p). The 
plotted errorbars also indicate the associated standard 
deviation (a^ r = y(4 — 7r)/(47rp)). When p reaches 
p S at, in few cases, the distance between two neighbour- 
ing nucleation centers is of the order of the nucleation 
distance. As a result, two centers can nucleate almost at 
the same time (avalanche effect) . The number of micro- 
cracks involved in such avalanches, also plotted in the 
figure, increases in the vicinity of p sat . This strongly 
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Fig. 18 Determination of the local propagation velocity c m 
of individual microcracks. Left: Evolution of the mean crack 
front position as a function of c m x t. Different curves cor- 
respond to different fractographic images with different mi- 
crocrack density. Slope of these curves defines the ratio A 
between continuum-level scale failure velocity v and microc- 
rack velocity. Right: Variation of c m with p. Horizontal red 
line indicates the mean value c m w 217m/s. 



suggests that the observed saturation in the d n evolu- 
tion (Fig. |l7) Top) results from this steric effect. 



5.3 Velocity of microcrack growth 

The fractographic analysis performed in Sec. [4] showed 
that, within the DZ, all microcracks grow at the same 
velocity c m . Unfortunately, it does not allow us to di- 
rectly measure the value of c m and its possible depen- 
dence with p. Indeed, Eq. [2] shows that the form taken 
by the conical brands depends only on the position of 
the nucleation sites and on the nucleation distance but 
not on the value of c m . 

We propose to use the deterministic reconstructions 
obtained in Sec. [4] to connect the macroscopic crack 
speed v to c m . Figure [18] A shows the evolution of the 
mean front position measured from reconstructions of 
different fracture surfaces (i.e. for different densities of 




2 4 

[Kx{l-v/C R )/y/r- 



■v/C d ] 2 (MPa 2 .m) 



Fig. 19 Variation of the fracture energy r with Kd = K x 
(1 — v/Cr)/ tJI — v/Cd- The vertical dashed line corresponds 
to microcracking onset K°j. = 1.83 MPay'm. Note that Kd can 
stricto-senso be identified with the dynamic SIF in absence 
of microcracking, i.e. for Kd < -RT5. The solid red line is a fit 
with Eq. [5] 



nucleation centers) as a function of numerical time step 
normalized by c m (homogeneous to a distance, corre- 
sponding to the distance travelled by a crack front at 
a velocity c m in absence of micracking events). All the 
curves are linear with a slope A that defines the ra- 
tio v/c m . At this point, one should recall that the re- 
construction procedure starts with an unrealistic initial 
straight crack front. We checked that this does not af- 



fect the measured value A (Guerra et al. 20121. Then, 
since v is known, c m can be deduced and plotted, on 
Fig. 18 B as a function of p. This microscopic velocity 
is found to be a constant, c m ps 217 m/s = 0.24Cr, 
independent of the loading K. the density p, and the 
macroscale velocity v. The implications of this result 
will be further discussed in the next section. 

The fracture energy measurements r(v) made in 
Sec. [3] will now allow us to understand why such a con- 
stant velocity c m is selected within the microcracking 
regime. Figure 19 shows the evolution of r with the 
dynamic SIF Kd- In the low speed regime (i.e. for low 
Kd) where crack propagates without involving any mi- 
crodamaging or microbranching, r scales with K\. We 
can also expect that the FPZ size R c scales with K\: 
R c = K^/auy where a is a constant of order 1, and ay 
is the material's yield stress. The scaling of F with Kd 
thus indicates a linear variation between -T and R r and, 



following Scheibert et al (2010), a constant dissipated 



energy e per volume unit within the FPZ. The fracture 
energy then writes: 



r 



— K l 



K. 



(5) 



Damage mechanisms in the dynamic fracture of nominally brittle polymers 



15 



As long as v < v a (thus without microdamage), one 
can relat e K d to K ( |Freund[ |1990[ ): K d = K x (1 - 
v/Cr)/ yj\ — v/Co where Cn refers to the speed of di- 
lational waves (cf. Tab. 1). By replacing K with K d in 
Eq. [T] and by subsequently injecting the result in Eq. 
[5j one gets the following expression for r vs. v. 



r 







1- 


eE 1-v/Cr p, 


ao\ 1 — v/Co 



(6) 



This expression yields a divergence of r at a finite ve- 
locity Coo given by: 



cE 



aa? 



- 1 



C R C 



D 



Cd — Gh 



(7) 



In this expression, the only unknown quantity is e/aay ■ 
This can be evaluated by fitting the first linear part (i.e 
for K d < K%) of the curve r vs. K^/E plotted in Fig. 
[191 with E q- 0ne then S ets: 



204 m/s = 0.23C R 



(8) 



This value sets the maximum crack growth velocity 
in the absence of microcracking: Beyond c^ dissipation 
diverges within the FPZ. Assuming that the same dissi- 
pation mechanisms are involved within the FPZ of each 
microcrack when v > v a , we expect that the limiting 
speed of microcracks is also determined by c^ , so that 



6 Role of microdamaging in the selection of 
continuum-level scale failure velocity 

The previous section has allowed us to better under- 
stand how the different variables that characterize the 
microdamage developing within the DZ are selected 
when v exceeds v a . The most surprising observation 
concerns the propagation velocity of microcracks c m 
that remains constant, of the order of 217 m/s, and 
significantly lower than the macroscopic velocity v. In 
other terms, damage does not slow down the macro- 
scopic crack tip as it was commonly believed until now 



(Ravi-Chandar 2004 Ravi-Chandar and Yang 1997 



Washabaugh and Knauss 19941, but on the contrary 



it accelerates it. And this acceleration factor is all the 
more important as p (or equivalently K) increases (Fig. 
20). 

The temporal evolution of a given point on the main 
front, shown in Fig. |2l] permits to better understand 




Nucleation center density (mm" ) 



Fig. 20 Acceleration factor A = v/c m as a function of the 
microcrack density p. Red line is a fit by Eq.[9]with b = 1.19± 
0.02. Horizontal dot line defines vt/c m and passes through 
the slope breakdown observed when p = p sa t (Taken from 
|Guerra et~al| ( |2012| ) 




Fig. 21 Position of a given point of the front as a function 
of "time" c m x t as obtained from the reconstructed simula- 
tion, applied to a fracture surface with p = 64.6 mm -2 . The 
dotted straight line indicates the expected position without 
microdamage (slope 1). The jumps (arrows) correspond to co- 
alescence events with microcracks (Taken from |Guerra et al| 
pH2| )) 



the effect of microdamage on the dynamics of macro- 
scopic cracking. The trajectory of this point exhibits 
random and sudden jumps. These jumps actually cor- 
respond to the coalescence of a microcrack with the 
main front. Between these events, the local front veloc- 
ity remains close to c m (i.e. the slope is equal to 1). 
Such an acceleration effect can be captured by a sim- 



ple mean-field model, detailed in Guerra et al (2012 1, 



which simply consists in counting the number of coales- 
cence events as the main front propagates over a unit 
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of length when the nucleation sites are located at the 
nodes of a square lattice. One then gets: 



.4 



1 



1 - bd n ^fp 



(9) 



where b is a numerical factor close to 1. This equa- 
tion, combined with the evolution of d n with p (red line 
in Fig. [17] Top) is plotted in red in Fig. 20 It reproduces 
quite well both the data coming from experiments (grey 
triangles) and from the reconstruction (black circle). 
Note the change in slope observed at p sa t- At this point, 
A = At, = 1.63. This corresponds to a ratio At, ~ Vb/c m 
where Vb — OACr is associated with the microbranch- 
ing onset. This strongly suggests that the steric effect 
at the origin of the saturation of d n also plays a role 
(still unsolved) in the micro-branching instability. 



7 Conclusion 

The series of experiments reported here were designed 
to shed light on the dissipation mechanisms that de- 
velop during fast crack growth in PMMA. They reveal 
that, above a well-defined critical velocity, PMMA stops 
to behave stricto senso as a nominally brittle mate- 
rial and crack propagation goes hand in hand with mi- 
crocracks that nucleate and grow ahead of the main 
front. Those let characteristic marks (referred to as 
conic marks) on post-mortem fracture surfaces, the mor- 
phological analysis of which allowed us to reconstruct 
the full dynamics of microfailure events. The simulta- 
neous space and time accuracies were of the order of 
the micrometer and the tens on nanoseconds, i.e. far 
better than what was reachable up to now. 

Analysis of the reconstructions demonstrated that 
the true local propagation velocity of single cracks re- 
mains limited to a fairly low value c m ~ 0.23cr in 
PMMA, while the apparent fracture velocity measured 
at the continuum-level scale (e.g. via potential drop 
method) can be much higher (up to twice larger!). Such 
an anomalously high measured velocity results in fact 
from the coalescence of microcracks with each other 
and with the main front - all of them growing at c m . 
In other words, the main effect of microdamaging is 
not, as commonly believed up to now ( |Washabaugh and] 
Knauss[ 1994; Ravi-Chandar 2004 Ravi-Cha ndar and 
Yang 19971, to slow down fracture by increasing dissi- 



pation within the DZ, but on the contrary to boost the 
propagation of the main crack. 

The value c m that limits the true local velocity of 
single crack fronts in PMMA in the dynamic regime is 
set by the physico-chemical dissipative and non-linear 
processes that develop within the FPZ, e.g. thermal 



(Estevez et al 2000), viscoelastic effects (Boudet et al 



1996[ Persson and Brener| |2005[ | or hyperelastic pro- 
cesses ( BouchbmdiretaIl|2008J|Llvne 2010 Buehkr 



et al 2003). Hence, it will depend on the considered 



material. On the other hand, the boost effect induced 
by microdamaging takes the form of a purely geomet- 
rical factor controlled by two quantities: the density p 
of nucleation sites and the mean distance d n at nucle- 
ation. These two variables fully characterize the damag- 
ing state within DZ and evolve with K. Once ascribed, 
they permit to relate the continuum-level fracture speed 
v to the true local propagation speed c m . Further work 
is required to unravel how PMMA internal structure se- 
lects the material-dependent quantities c m , p{K), and 
d„. 

The presence of microcracks forming ahead of a dy- 
namically growing crack has been evidenced in a vari- 



ety of materials, e.g. in most brittle polymers (Ravi- 



Chandar 1998 Du et al 20101, in rocks (Ahrens and 



Rubin 1993 1, in some nanophase ceramics and nanocom- 



posites (Rountree et al 20021, in oxide glasses (Roun- 



tree et all 2010 1, in metallic glasses (Murali et al 



2011), 



etc. The boost mechanism demonstrated here on PMMA 
is expected to hold in this whole class of materials. Note 



that a similar effect has also been invoked (Prades et al 



2005) in oxide glasses under stress corrosion, where 



ultra-slow crack propagation (down to few tenths of 
nanometers per second) was found to involves "nano"- 
cracks forming ahead of the main crack tip ( |Celarie| 
et~aH [20031 |Celarie et al[ [2003] |Bonamy et al[ |2006| 
Ferretti et al 20111. This yields us to conjecture that 



failures with anomalously high apparent velocities mea- 
sured at continuum-level scale may arise in all situa- 
tions involving propagation-triggered microcracks, in- 



cluding e.g. shear fracture in compressed granite (Moore 



and Lockner 1995 1 , thermal failure in shale (Kobchenko 



et al 2011 Panahi et al 20121, and more generally 



failure of so-called quasi-brittle materials. One of the 
most robust observations in this field is the power-law 
form followed by the distribution in size and in time 
interval between two successive microfailure events (see 
(Bonamy 2009 Deschanel et al 20091) for recent re- 



views) . It would be interesting to see whether or not the 
microcracking evidenced here in dynamic failure regime 
share the same scale-free features and how this affects 
the boost effect. Work in this direction is under way. 
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